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Abstract. In this paper we introduce a new class of integrable systems, naturally associated to 
Hurwitz spaces (spaces of meromorphic functions over Riemann surfaces). The critical values of the 
meromorphic functions play the role of “times”. Our systems give a natural generalization of the 
Ernst equation; in genus zero they realize the scheme of deformation of integrable systems proposed 
by Burtsev, Mikhailov and Zakharov. We show that any solution of these systems in rank 1 defines 
a flat diagonal metric (Darboux-Egoroff metric) together with a class of corresponding systems of 
hydrodynamic type and their solutions. 

1 Introduction 


Deformations of Riemann surfaces appear in different aspects of the theory of integrable systems. 
We mention the theory of systems of hydrodynamic type (in particular, the theory of Whitham 
deformations of integrable systems) nziia, theory of algebro-geometric solutions of equations with 
variable spectral parameter [Tlj (the main representative of this class is the Ernst equation from 
general relativity) and the theory of Frobenius manifolds pm. 

The Ernst equation 

((? - + ((^ - OG^-G-i)^ = 0 , 

where G{^, is a matrix-valued function (the stationary axially symmetric Einstein equations are 
equivalent to this equation if G is a 2 x 2 matrix with some additional symmetries) is the compatibility 
condition of the following linear system [U — V pair) 122]: 


G.G-^ ^ 9T GfG-i ^ 

— = —5-T , ^ ^-T 

1 — u 1 + V 


( 1 . 1 ) 


where v is the following function of spectral parameter I G C and variables (CjO- 
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Function iy{l) is nothing but the uniformization map of the genus zero Riemann surface which is the 
twofold covering of Z-plane with two branch points at / = ^ and I = The map l{i') is a rational map 
—> CP^ of degree two with critical values ^ and 

If the matrix dimension of G is 1, we can introduce the function / = In G and the Ernst equation 
turns into the Euler-Darboux equation: 


d^f I df 1 df 


( 1 . 2 ) 


The natural question is whether the Ernst equation is an isolated example of an integrable system 
related to the space of rational maps or it is possible to define natural analogs of the Ernst equation 
which would correspond to spaces of rational maps of arbitrary degree? More general, is it possible to 
go beyond the spaces of rational maps, and define natural analogs of the Ernst equation corresponding 
to general Hurwitz spaces of meromorphic functions of degree N on Riemann surfaces of genus 
gl Is there any link between these higher analogs of Ernst equation and existing theories of systems 
of hydrodynamic type mi El, Frobenius manifolds and Darboux-Egoroff metrics corresponding to 
Hurwitz spaces 10]? 

If we assume that v in in]) is a constant, then the compatibility conditions of mi leads to 
equation of principal chiral model; therefore, according to terminology proposed by Burtsev, Mikhailov 
and Zakharov E], it is natural to call the Ernst equation the “deformation” of the principal chiral 
model equation. In 10] it was studied the general problem of deformation of a given integrable system 
which has U — V pair, where matrices U and V are meromorphic functions of constant spectral 
parameter ly. If one allows ly to depend on space variables {x,y), then the zero curvature condition 
Uy — Vx + [U,V] = 0 implies a set of differential equations for “variable spectral parameter” u and 
poles of matrices U and V. However, in 10] no regular method was given to solve these differential 
equations. 

The first goal of this paper is to fill this gap. Namely, we show that the deformation scheme of 
10] can be realized in terms of the spaces of rational maps of any given degree N. In this way 
we get a new hierarchy of non-autonomous nonlinear integrable systems. We show how solutions of 
the new systems can be described in terms of matrix Riemann-Hilbert problem and isomonodromic 
deformations. If the matrix dimension equals 1, these non-linear systems give rise to systems of linear 
non-autonomous second order partial differential equations, generalizing the Euler-Darboux equation. 

Second, we extend our framework to construct a class of new integrable systems starting from an 
arbitrary Hurwitz space (space of meromorphic functions of degree N on Riemann surfaces of 

genus g) for g >2. 

Third, in rank 1 (by rank of the integrable system we mean its matrix dimension) and any genus 
we observe a very close relationship between our systems, Darboux-Egoroff metrics, systems of hy¬ 
drodynamic type and Frobenius manifolds. Moreover, our general formalism allows to give a simple 
description of systems of hydrodynamic type (as well as their solutions) associated to Hurwitz spaces 
Hg,N- 

Let us describe our results in more details. Each rational map R{'j) of degree N defines an Wfold 
branch covering L of CP^; a point P of L is a pair (^, 7 ) such that 


I = R{i) • 


(1.3) 


We assume that R{’j) maps the infinity in 7 -plane to the infinity in /-plane and P( 7 ) = I + o(l) as 
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7 —> oo; then function Rij) has the form 


N-l 

Rii) = 'y+Y^ 

m=l 


^ m 
1 - 


(1.4) 


Denote the critical points of the map i? by 71 ,... ,7m; we shall assume that all of them are simple; 
then M = 2N — 2 according to the Riemann-Hurwitz formula. 

The ramification points of the covering L are denoted by = (Irmlm) (tliey are simple since oil 
critical points of R{'j) are simple); their projections Im on ?-plane are called the branch points of the 
covering L (we adopt terminology of i); these are the values of the rational map R{^) at its critical 
points: Im = Rijm)- In the sequel we shall assume that all Im are different. 

Introduce two functions on the covering L: the function vr, which projects L on /-plane: vr(P) = /, 
and function z/, which projects L to 7 -plane: ^{P) = 7 . The map z/ : L —> CP^ is a one-to-one map; 
its inverse is nothing but uniformization map of the covering L. The maps i' and vr are related as 
follows: R{i'{P)) = ■ 7 r(P). 

Due to our assumption about the behaviour of ^( 7 ) at infinity, in a neighbourhood of the infinite 
point on some (we shall call it the first) sheet of L the map iy{P) behaves as follows: v{P) = / -|- o(l). 

The structure of Riemann surface on the branch covering L is defined as follows: in a neighbourhood 
of any point where L is non-ramified we can consider I as local parameter. In a neighbourhood of a 
ramification point P^ the local coordinate is chosen to be Xm = \// — Im- 

The branch covering L is completely defined by the branch points Im and a representation of 
fundamental group 7 ri(CP^\{/i,... ,/m}) in permutation group Sn- An element of permutation group 
describes permutation of sheets of the covering L if / encircles a given contour in CP^ \ {/i,... ,/m}- 
We shall consider local deformations of the covering L such that this representation is kept fixed. Then 
the branch points Im can be considered as natural local coordinates on the space of rational maps; 
they will play the role of independent variables of our systems. 

Let us fix some point Pq of L such that its projection on /-plane /q = vr(Po) is independent of 
all {/m}; let 7 o = i^(Po)- Consider the following system of matrix linear differential equations for an 
auxiliary matrix-valued r x r function 'L(/,{/m}): 


dim 


7o - 7m 

^{P) - Im 




We assume that the solution 'L of (USD is normalized by the condition 


(1.5) 


$(/ = oo) = /. (1.6) 

If we put in (ESI) P = Po, then the compatibility conditions of the system (ESI) imply that all Jm 
are logarithmic derivatives of a matrix-valued function G{{lm})'- Jm = Besides that, the 

compatibility conditions of the system ESI imply the following coupled system of non-autonomous 
non-linear matrix partial differential equations of second order: 

((70 - = ({70 - 7..)0,.G-'),_^ (1.7) 

for all m, n = 1,..., M. We call the matrix dimension r of the matrix G the rank of the system E2I). 

For rational maps of degree two the system E21) coincides with the ordinary Ernst equation EU). 
If the rank equals 1, the systems EH) gives rise to certain generalization of the Euler-Darboux equation 

E3. 
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We define the tau-function of these systems by the following system of compatible equations 


A. 

dlrr, 


Inr = —res 
2 


tr 


I {dp^^ 


- 1\2 


dl 


( 1 . 8 ) 


where dp^ = {d^/di'{P))di'{P); these equations fix r up to an arbitrary constant multiplier if we 
assume that r is holomorphic function of {Im}- 

Each system GUI possesses a subclass of “isomonodromic” solutions which can be built from an 
arbitrary solution of the Schlesinger system in the same matrix dimension. For this isomonodromic 
sector of solutions we link the tau-function to the Jimbo-Miwa tau-function of the Schlesinger system. 
This link generalizes the formula (found in |15j i relating the so-called “conformal factor” of the Ernst 
equation with Jimbo-Miwa tau-function. We also show how to construct solutions of the systems jni) 
from solutions of matrix Riemann-Hilbert problems on CP^. 

After developing the theory of the systems cm, which are related to the space of rational maps 
Hq^n, we formulate a similar construction based on an arbitrary Hurwitz space 

Consider in detail the rank 1 case. Let the A^-fold branched covering L have genus g; then the 
number of branch points is equal to M = 2g + 2N — 2 (as before, we assume that all the branch points 
are simple and have different projections on Z-plane). Again, the projections Im of the branch points 
Pm on /-plane can be used as local coordinates on Hg^pf. Let us introduce some basis of canonical 
cycles on L. Denote by E{P,Q) the prime-form on L, and by Bm the normalized (all a-periods 
vanish) abelian differential of the second kind with the unique pole of second order at Pm with leading 
coefficient 1. Let us choose two arbitrary points Pq and Qq on L such that their projections 7r(Po) 
and 7t(Qo) on /-plane are {/m}-independent. Then the Bg,N analog of the genus zero systems (II.7|) in 
rank 1 looks as follows: 

b{Pm,Pn) f Vn df ^ Vm Of 

din 


dlm.dln 


dl 


H- 

m 


= 0 , 


m ^ n 


(1.9) 


where 


b{Pm,Pn) = 


BiP,Q) 


P = Pm,Q = Pn 


dXm(,P^^^n{Q^ 

and B{P, Q) = dpdq lnE(P, Q) is the Begmann kernel; Vm = m • 

Solutions of the system dn can be constructed as follows. Let / be an arbitrary closed contour 
on L such that its projection on the /-plane 7r(/) is independent of {Im} and Pm 0 / for any m. Let 
h{P) be an arbitrary independent of {Im} Holder-continuous function on /. Then the function 

E{Po,P) 


f = (ph{P)dp In ■ 

Jl 


( 1 . 10 ) 


E{Qo,P) 

satisfies the system dm. 

The systems GSD turn out to be in close relationship with diagonal flat metrics (Darboux-Egoroff 
metrics). Namely, consider the diagonal metric 


M 

ds^ = 

m=l 


Qmmdl, 


2 

m 5 


(1.11) 


where gmm = dlnr/dlm and r({/m}) is the tau-function corresponding to an arbitrary solution of 
dISI). The rotation coefficients of this metric turn out to be given by the Bergmann kernel: 


Pn 


= 7;b{Pm,Pn) 
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i.e. they depend only on the covering L and don’t depend on a particular solution of (jl.DI) . In 
accordance with the Rauch variational formulas |H], which describe dependence of the Bergmann 
kernel on the branch points, these coefficients satisfy the equations 


dh 




n 


for distinct I, m, n, which, together with annihilation of all (3mn by the operator 9/dlk , guarantee 

flatness of the metric cm- 

The Darboux-Egoroff metrics are known to be closely related to integrable classes of systems of 
hydrodynamic type, including, in particular, Whitham equations which describe slow deformations of 
the Riemann surfaces arising in dispersionless limit of finite-gap solutions of integrable systems. These 
systems are solvable via so-called generalized hodograph method EHIIIII- Our present scheme gives 
a short and simple formulation of the theory of the systems of hydrodynamic type corresponding to 
Hurwitz spaces. 

We notice that the present paper has some overlap with the paper d, devoted to isomonodromy 
deformations in higher genus, which appeared simultaneously with the first version of this text d- 

The paper is organized as follows. In section 2 we derive an auxiliary system of differential 
equations which describe dependence of the critical points of the rational map of the form JLI on 
the critical values. Then we show integrability of the systems (El), define corresponding tau-function, 
and discuss the relationship of these systems to the Riemann-Hilbert problem and the Schlesinger 
system. Here we also show how the systems (El are related to deformations of integrable systems 
proposed in jH]. In section 3 we define analogs of the systems (ll.7|l related to Hurwitz spaces Hg^j^ 
in arbitrary genus g > 2, and discuss their properties. In Sect.4 we show that each solution of these 
systems in rank 1 defines a Darboux-Egoroff metrics, and discuss related systems of hydrodynamic 
type in our framework, together with their solutions. Here we also outline the link between higher 
genus analogs of systems (El and isomonodromic deformations on algebraic curves. In section 5 we 
discuss potential directions of future work. 


2 Non-autonomous integrable systems related to spaces of rational 

maps 

2.1 Differential equations for critical points of rational maps 

Consider a rational map R{'y) of degree N of the form (|1.41) . If we choose the critical values li,... ,1 m 
of the map (El as independent parameters, each critical point 7 ^ becomes a function of all {Im}] 
the study of these functions is the main subject of this subsection. The function ^{P) depends on the 
variables li,... ,1m as parameters. In the sequel we shall denote the point of L which belongs to the 
jth sheet and has a projection I on the Z-plane by The map ^{P) has its only pole at the point 
at infinity of some sheet of L; we enumerate the sheets of L in such a way that this sheet has number 
one; therefore, 

v{P) = l + o{l) as P ^ (2.1) 

Theorem 1 Function n{P), considered locally as function of I and depending on branch points li,... ,1m 
as parameters, satisfies the following equations: 

^ an dv _ an 

dl ^ - 7 n ’ dl, 


5 


i^-7n 


( 2 . 2 ) 






where a^, are some functions of the branch points. 


Proof. Consider the local behavior of the function y{P) near a branch point: 

— Tm “1“ "x/^ l-m P 0(1 Im) S-S P ^ Pm ■ (^-3) 


From (lOl) we conclude that 


du 

'm 


^m 

2^1 — Ijn 


+ 0{l) , 




+ 0 ( 1 ) , 


as P —> Pm- By H2.dl) . we can rewrite these expansions using z^(P) as the global coordinate on L: 


dn 

'm 


2(z^ Tm) 


+ 0 ( 1 ), 


dn 

din 


-dn 


'2(l^ - 7n) 


+ 0 ( 1 ) 


(2.4) 


where 7 ^ = v{Pn)- Moreover, from H 2 . 1 B we conclude that dv/dl = 1 + o(l) and du/dln = o(l) as 
n{P) — > 00 . Therefore, dv/dl is a meromorphic function of n with simple poles at all the points 7 ^ 
with residues k ^/2 and value 1 at infinity. Analogously, the function dv/dln is a meromorphic function 
on CP^ with simple pole at In and zero at infinity. Therefore, we get equations H2.2|) with an = k^/2. 
□ 


Corollary 1 The critical points 'jm of the rational function and residues am from ifOl) depend 
as follows on the critical values Im-' 


d^m _ 

din 'Jn 7m 


dam ‘2-OinO^m 
din (7n - 7m)^ ’ 

for all m,n = 1,..., M. 


m n ] 


dim 

dim 


M 

1+ E 

n=l,n^Tn 


CXri 

7m - 7n 


m n \ 


dam 

dim 


M 

E 

72 = 1 , 


{'In ~ 7m)^ 


(2.5) 


( 2 . 6 ) 


Proof. Equations (EH) and ESI follow from compatibility of equations (I2.2|l . 

□ 

Rational functions of the form EH were introduced by Kupershmidt and Manin m in connection 
with Benney systems. The fact that the critical points of these functions satisfy equations (IT31).(I^. 
follows from the recent paper by Gibbons and Tsarev m However, we did not find in existing 
literature the whole set of equations EH) (12.51) . (|2.6I) associated to the functions EH- 

Two-fold coverings. For an illustration consider the case N = 2, when the covering L has two 
sheets and two branch points li, I 2 . Then the rational function P( 7 ) (II .41) can be explicitly written in 
terms of its critical values li and I 2 '. 


Ril) = 7 + 


(^1 — h) 
16(7 - H^) 


^We thank E.Ferapontov, who attracted our attention to this work 


(2.7) 
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The map ^{P) looks as follows: 


v{P) “ 2 ^ 2 a/(^ “ h){l - h)^ ■ 


The critical points 71^2 and variables ai ^2 are given by: 


71 = ^{h) 


3li + I 2 
4 


72 = u{l2) 


ll + 3/2 
4 


( 2 . 8 ) 


(2.9) 


_h-h 

0^1 — — 0-2 — - 


( 2 . 10 ) 


2.2 Spaces of rational maps and new hierarchy of non-antonomous integrable 
systems 

Starting from an arbitrary branch A^-fold covering L of genus zero we can construct a hierarchy of 
integrable systems as follows. 

Fix some point Pq £ L such that its projection Iq on CP^ is independent of all {Im}, be. 70 = i^{Po) 
depends on {Im} according to the equation 

P{7o{h, ■ ■ ■, Im)) = ^0 ■ 


Consider the following system of first order differential equations for a r x r matrix-valued function 

^{PAlm})-- 


_ 7o - 7m 
dim v{P) - Jm 


( 2 . 11 ) 


where Jm are r x r matrix-valued functions of {Im}- 

As a corollary of compatibility conditions of the linear system dTTTl) functions Jm can be expressed 
in terms of the single function G = 'I'(Po): 


^ra — 07 

^I'm, 


( 2 . 12 ) 


moreover, the function G satisfies the following system of non-linear partial differential equations: 


iilo-lm)Gi^G-^)^ ={{jo-ln)Gi„G-^)^ . (2.13) 

(to derive (12.1311 from (tTTTl) one needs to make use of equations (EH, (ESI) and (EH). 

We write the “total” derivative of T with respect to Im in Trm to emphasize that in (CTl) we 
consider I and li,. ■ ■ ,lm as independent variables. 

Alternatively, consider 'k(T’) as a function of ^{P) and Ii,.--,Im- Then, since v{P) is itself a 
function of I and {Im}, we can rewrite (j2.1 1 11 using the chain rule 

dT _ 5T du aT 
dim dim dim du 

where the notation /dim means that z^(P) remains fixed. 

Then, using (I2.2|l for du/dlm, we rewrite (j2.11ll as follows: 


(9T(P) am d^{P) _ -fo-jm j 
dim u{P)-Jmdiy{P) u{P)-Jm " ^ ^ 


(2.14) 
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A simple calculation making use of equations and system (izni), dsni), shows that if G{{lm}) 

is a solution of non-linear system (EISl), then the 1-form 


M 


q = 


(70 - 7m)^^ 

^Lxr, 


m=l 




(2.15) 


is closed, dq = 0. 

The closedness of the 1-form q implies the existence of the potential r , which can naturally be 
called the tau-function of the non-linear system 


Definition 1 The function r(/i,..., Im), defined by the following system of equations: 

51nr _ (7o-7m)\„^^ ^-n 2 
—2^— ) 

uim 

up to an arbitrary constant multiplier, is called the tau-function of integrable system JiKW . 


(2.16) 


Using equations (irni) o, Q, we find that the second derivatives of the tau-function are given 
by the following expression: 


S^lnr _ ( 70 -7m)(7o-7n),,, 

Ml” 2h„-7„)2 > 


(2.17) 


for m ^ n. 

Taking the residue of the linear system at P = Pjn, we find 

= (70 - 7m)Gi^G-^ ; 

due to this relation the definition allows an alternative formulation: 

Definition 1’ The tau-function of the system JE71) is defined by the following equations: 


(2.18) 


Slnr 1 


dir. 


where 






I tr(dpT 'k 


- 1\2 


dl 


-1, 


(2.19) 


( 2 . 20 ) 


Let us prove the equivalence of the two definitions of the tau-function. Using ^{P) as a global 
coordinate on L, we have: 

dl = T—di' ; 
ah' 

therefore, 


tr{(dpT^-i)^} du 


M 


dl 


= —tr{(T,M/-i)2}dj. = <1 1 + ^ 


Or 


tr{(^,^'-i)2}dj. , 


m=l 


^ -Im 


and (ITT^ coincides with (HIHI). 












The Ernst equation. For N = 2 the hierarchy (l2.i;-iB reduces to a single equation. If one chooses 
point Pq to coincide with [i.e. the point of L where I = oo and in a neighbourhood of which 

= -1 + + / 2)/2 + o(l)], we get 


h + h 

70 = u{Po) = . 

Taking into account expressions (EHl), we have 

h — h h — h 

70 - 71 = —^— , 70 - 72 = —^— • 

If we now assume that h and I 2 are conjugated to each other: li = I 2 
takes the following form: 


( 2 . 21 ) 

then equation 


((e - + ((^ - = 0 . (2.22) 

This equation is called the Ernst equation; it is equivalent to vacuum Einstein’s equation for stationary 
axially symmetric spacetimes (the matrix G in this case must be real, symmetric and must have unit 
determinant). In this case the linear system (12.111) is equivalent to the Lax representation of the Ernst 
equation found in 1978 by Belinskii-Zakharov [2j and Maison m- We notice that the Maison’s Lax 
pair is a partial case of our linear system written in the form (EUl); whereas the Belinskii-Zakharov 
linear system is a partial case of (EH. 

Due to expressions for ai^ 2 , the definition of the tau-function r can be written down as 

follows: 

( 0 , 0-7 . ^ . ( 2 . 23 ) 

Formula coincides with the definition of the so-called conformal factor - one of the metric 

coefficients which correspond to a given solution of the Ernst equation. 


2.3 New integrable systems and deformation scheme of Burtzev-Mikhailov-Zakharov 

The possibility to construct the class of “deformed” integrable systems, or integrable systems with 
“variable spectral parameter”, different from the Ernst equation, was first discovered by Burtzev, 
Mikhailov and Zakharov |3j. They proposed to consider the Lax pairs of the form 

^ = ^ , (2.24) 


where x and y are independent variables; matrices U and V depend on (x, y) and the “variable spectral 
parameter” v (which in turn depends on (x,y) and the “hidden” spectral parameter /): 


Ni 


TT/ \ / \ Un{x,y) 


N 2 


^ \ \ Vn{x,y) 

V{x, y, v) = Vq{x, y) + > —r 

^ z^(x,y) -7„(x,y) 


(2.25) 

(2.26) 
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As a part of compatibility conditions of the linear system (IT^ . after an appropriate fractional-linear 
transformation in z^-plane, the following system of equations for u{x, y, 1) must be satisfied: 


du 

dy 


N2 

+ E 

m=l 


bm 

1^-7m 


= 0 , 


dv 

dx 


+ E 


m=l 


Cm 

I' -7m 


= 0 , 


(2.27) 


where bn and Cn are certain functions of (x,y). The compatibility condition of the system H2.27I) gives 
the following system for 7 „(x,y) and 7 n{x,y)\ 


d7n 

dy 


-f 


N2 

E 

m=l 


In '7m 


= 0 , 


9% 

dx 


Ni 

“ 7n - 7m 

m=l 


= 0 , 


(2.28) 


f)c 

uon 

dy 


— 2Cr 


N2 

E 


(7n 7m)‘ 


= 0 


dbn 

dx 


Ni 


-2bnY, 


(7n 7m)^ 


= 0 


(2.29) 


It is easy to establish the relationship between solutions of the system (ITTtI) . (IT^ . TTM . and 
solutions of our system (lO) . (IT^ . (EH). 

Namely, suppose that function {lm}m=i) satisfies equations EH v^ith respect to variables /m* 
Assume that M = A'i-I-A '2 and split the set of variables {/i,..., /Ari+Ar 2 } into two subsets: {li,... ,1^^} 
and {fi,..., where In = Ini+u, n = 1,..., A^ 2 - In the same way we split the set { 7 m} of values of 
function I'iP) at these points: 


{7i,---,7m} = {7i,...,7Ari}U{7i,...,7Ar2} , 


where % = 7Ni+n, n = 1,..., N' 2 - 

Now assume that the “untilded” variables are arbitrary functions of variable x and 

the “tilded” variables are arbitrary functions of variable y. Then using (EH we get the 

derivative of ^{P) with respect to x: 


9v _ 9v dim 
dx ^ dim dx 

m=l 


E 


m=l 


dir. 


OLr 


dx 


(2.30) 


therefore, 


dv 

dx 


Ni 

+ E 

m=l 


Cm 

y -Im 


= 0 , 


(2.31) 


where Cm = this coincides with the second equation in (I2.27|) . The first equation in (I2.27|) is 

obtained in the same way after identification 


bm — OiNi+m 


dlNi+m 

dy 


Equations (ItM and (IT^ for 7 „, bn and Cn as functions of (x, y) arise as compatibility conditions 
of equations for and Vy. 

Therefore, spaces of rational maps of given degree provide solutions of the system (iT^Tll . (ItM . 
TTMi if 


1. We split the set of the branch points {Im} into two subsets and 


2. Assume that one subset contains the branch points which are arbitrary functions of x only and 
another subset contains the branch points which are arbitrary functions of y only. 

Concluding, we see that the system (I2.13|) provides a realizations of the deformation scheme of 0. 
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2.4 Solutions via the Riemann-Hilbert problem 

Solutions of systems (TOl) in terms of matrix Riemann-Hilbert problem are given by the following 
theorem. 


Theorem 2 Consider a closed contour I on the branch covering L such that its projection tt{ 1) on the 
l-plane is independent of {Im}- Assume that none of the ramification points Pm belongs to 1. Define a 
non-degenerate matrix function I 3 P ^ H{P) which is independent of {Im}- Suppose that a function 
4'(P) satisfies the following Riemann-Hilbert problem on L: 

1. 4'(P) is holomorphic and non-degenerate on L outside of the contour I where it has the finite 
boundary values related as follows: 


4'+(P) = ^_{P)H{P) . 


(2.32) 


2. 4* satisfies the normalization condition at (this is the point of L such that 7r(oo^^^) = oo 

and = oo): 

4'(oo(^)) = I. (2.33) 

Then the function 'h satisfies the linear system i2.11l) with 

G = 4/(Po) , (2.34) 

and, therefore, the function G solves the system 

Proof. Let 4' satisfy the Riemann-Hilbert problem. Consider its logarithmic derivative 4'/^ 4^“^. 
Due to the independence of the contour I and the matrix H[P) of Im, this logarithmic derivative is 
single-valued on L. Moreover, it is obviously holomorphic on L outside of the point Pm- Let us write 
the first two terms of the Taylor expansion of 4'(R) near Pm'- 


4 - = 4-0 + \// - Im^i + ... ; 

then 

4/;^4/-i = + 0(1) . 

Zy L Ljyi 

Therefore, in terms of the uniformization map ^{P), we can write 


4'/ 4'"^ 

I'm 


G„ 


J^{P) - Im 


with some matrix Gm which is independent of P. The constant term in this formula is absent due to 
normalization condition (|2.33|) . To compute Cm we put P = Pq, i.e. ^{P) = 70 ; then, using (|2.34|1 . 
we get 

Cm = (70 — 7m)Gi^G ^ , 

which shows that the function 4^, indeed, satisfies (tTTTl) . and the corresponding function G solves the 
system 


Below we also establish a relationship between systems (ITT^ and another type of the Riemann- 
Hilbert problems, where the function 4' is allowed to have regular singularities. 
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2.5 Relationship to isomonodromic deformations 

The set of solutions of each system (imi has a subset of solutions corresponding to isomonodromic 
deformations of ordinary differential equations with meromorphic matrix coefficients: 

f = P-35) 

where Aj £ gl{r) are certain matrices which are independent of 7 and such that ~ Consider 

GL(r)-valued solution ^( 7 ) of (12.,‘1511 satisfying the initial condition at some point 70 £ CP^: 

T(7 o) = I . (2.36) 


The solution ^( 7 ) has regular singularities at the points Zj] this function is generically non-single- 
valued in the 7 -plane: it gains the right multipliers Mj under analytical continuation along contours 
starting at 70 and encircling poles Zj. The matrices Mj are called the monodromy matrices of equation 
((TSSJ); they generate the monodromy group of equation (ESSJ). If all the monodromy matrices are 
independent of the positions of singularities {zj}, then (in the generic case, when none of the eigen¬ 
values of each matrix Aj differ by an integer number) the function T satisfies the following equations 
with respect to the positions of singularities Zj: 




A. 


dzj V 7 o - Zj 


J-Zj 


T . 


(2.37) 


Compatibility of equations and (lOTll is equivalent to the Schlesinger system for the functions 

^j{{Zk})- 

dAj [dj, dfc] [Aj, Ajf] ■ / 1 

= — -^- , J 

dzk Zj - Zk 70 - Zk 


dzj 



[dfc, dj] \ 

Zk - 70 J 


(2.38) 


The tau-function tjm of the Schlesinger system, introduced by Jimbo, Miwa and their collaborators 
is defined by the following system: 


d 


— InTJM = -res|^=^.tr 


Qtjm 

dzj 


(2.39) 


Each solution of the Schlesinger system induces a solution of hierarchy (TOl) according to the 
following theorem: 


Theorem 3 Consider a solution {Aj{{zk})} of the Schlesinger system h2.St^) . together with the tau- 
function Tjm and the corresponding solution ^( 7 , {zj}) of the linear system )2.H5\) . \2.3T^ normalized 
by h2.36[) . Let L be a genus 0 covering of the l-plane with simple ramification points Pi,...,P m/ 
choose the map ^{P) to satisfy \2.1\) . Consider an arbitrary set of points Pq and Qi, ■ ■ ■ ,Ql on L 
such that their projections Iq and gi ,... ,/xt on CP^ (respectively) are independent of the branch points 
li,... ,1 m■ Let us assume that the arguments of the solution and the tau-function of the Schlesinger 
system are given by the formulas 


Zj = u{Qj) , 70 = i/(Po) . 


(2.40) 
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Then the function 


^ (7,7o,{^i}) 


-y=u{P),'yo=u{Po),Zj=u{Qj) 

satisfies the linear system nn of the hierarchy JE73I) with the functions Jm defined by 






7m -70 


■^( 7 ?) 


(2.41) 


(2.42) 


Therefore, according to mn), we have Jm = Gi^G ^ for some function G{{lm}), which satisfies the 
system 1123). 

The corresponding tau-function r is related to the Jimbo-Miwa tau-function as follows: 


i=i 


trA2/2 




TjM{{Zj}) 


Zj —i.'(Qj) 


(2.43) 


where the derivative dv/dl is given by equation hJJA) . 

Proof. Taking the derivative of the function 'k(P) with respect to Im using the chain rule, we get 

L 


5T(P) _ (9T du{P) ^ ^ ^ 


dir. 


dv{P) dir 


dzj dljyi 

J=1 3 




OLr, 






fZi ~ jZi 


-T 




E 


= 7. j„<i, 


7m 7m ^(-P) 7m 

where functions Jm are defined by (12311). 

Let us show how to prove the relation (IO 3 I 1 between tau-functions. Taking into account the 
definition of Jimbo-Miwa tau-functions (lOnil . we have: 

Now, using definition (ITTO of the tau-function r, we get 

L 

\P=Pr, ~ 


9 In r a, 


-it2 


dim 2 
By (12.21) . we see that 


= 7^tr(M-i) 


(Xr 


tr^2 

j^l ^ (7m “ Zj)"^ ^ 7m — 


+E 


dzi In TjM ■ 


(2.44) 


dl 




(Xr. 


dzj 




{jm-Zjfi' dim Im-Zj 


therefore, applying the chain rule in (12.441) (and taking into account that are integrals of the 
Schlesinger system), we come to (12.431) . □ 

Remark 1 The relationship between the systems (tniTl) and isomonodromic deformations is a gener¬ 
alization of the link between the Ernst equation and the Schlesinger system established in m 
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2.6 Rank 1 systems 

When the function G from is a scalar one, the system (im can be rewritten as a system of 

linear scalar second order differential equations in terms of the function = InG: 


(7m In) 


dlm.dlr 


CXn 


7m - 70 df 

7r 7o df'm. 




7n - 70 df 
7m - 70 din 


= 0 , 


m ^ n 


(2.45) 


In derivation of the system (inSD from (tTTlTl) we used equations (12.51) . 

In particular, any solution of the matrix system gives a solution of the scalar system 

if we put / = In detG. 

The linear system (CTl) turns in rank 1 into the scalar system 


dip{P) 7o - 7m df 


dir. 


u{P) - 7m dir 


(2.46) 


where '(/’(P, {^m}) = la 'h. As well as in the matrix case, the function ip is generically non-single-valued 
on L. 

The definition of tau-function (ItTgI) now looks as follows: 


51nr _ (70 - 7m)^ f df 

dl m. 2am. 1 dl m 


Alternatively, it can be rewritten in terms of the function 'ipiP), using (12.191) : 


(91nr 1 

—-— = -res 

dim. 2 


r 

P" \ dl 


(2.47) 


(2.48) 


where d'ip{P) = 'ip^{P)di'{P). 

Let us discuss the solutions of system (ITIHI) and equation dunj). 

Theorem 4 Let I be an arbitrary smooth closed contour on L sueh that its projection on the I-plane 
Tr{l) is independent of {Im} and Pm 0 I for any m. Consider on I an arbitrary Holder-continuous 
function h{P) independent of{lm}- Then the function 

/=/«« (2,49) 

Ji i^{Q) - 70 

satisfies system i2.45\ ). 

Proof. In the scalar case we know explicitly the solution of an arbitrary Riemann-Hilbert problem 
from theorem [21 If we define the function 

h{P) = ^. In H{P) 

(we assume that it is single-valued on I i.e. the index of H{P) on I equals 0), the Riemann-Hilbert 
problem firm becomes additive: 


f;+{P) = TP-{P) + 27rih{P) 


(2.50) 
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for P G I (we recall that V' = In'h); here V'i stands for the boundary values of the function ^ on 1. 
Normalization condition turns into 

^(oo(^)) = 0. (2.51) 


The solution of (lT50ll . (ITSTI) for a Holder-continuous functions h{P) is given by the Cauchy integral 


V>(P) = 


/ h{Q)dv{Q) 
v{Q) - v{P) ’ 


which implies at P = Pq, when z^(P) = 70 - n 


(2.52) 


Formula (inuii can also be verified as follows. First, one can check by the direct substitution, using 
equations (fO) . (IT31) . (fT^ . that the function 


1 dv{Q) 
^{Q) - 70 dl 


(2.53) 


satishes system (iTlSl) for any I = 7r{Q) independent of {/m}- Using linearity of equation (1051) we can 
consider the superposition of these solutions at different Q G I with an arbitrary {/m}-independent 
measure h{Q), which leads to (I2.49|) . 

In fact, we can consider any (say, compact) subset P C L whose projection on Z-plane is {Im}- 
independent, and such that Pm 0 D. We can define an arbitrary {/mj-independent measure d^{Q) on 
P; then the superposition principle implies that the function 


is a solution of ( 1 ^^ . 



h{Q) dv{Q) 
^{Q) - 70 dl 


d^Q) 


(2.54) 


The Euler-Darboux equation. As we noticed earlier, for the twofold covering with two branch 
points ^ and h = i the system is equivalent to Ernst equation (EHU). In scalar case we 

get the following equation in terms of / = In G: 


f--APi =0 

% 2 ({-{) 

If we introduce the real coordinates {z, p) such that ^ = z + ip, ^ = z — ip, this equation takes the 
form of Euler-Darboux equation: 

fzz H —fp + fpp = 0 . (2.55) 

P 

From ( 1 ^ . taking into account that Pq = 00 and 70 = (/i -|- we obtain: 


1 du{P) _ 1 

HP)-JO dl - [(/-^)(Z-e)]V 2 • 

Then representation gives the solution [3] of the Euler-Darboux equation: 


/ 



h{l)dl 
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3 Generalization to higher genus 

Here we discuss possible ways to define systems, analogous to (Eini), starting from Hurwitz spaces 
Hg i^ (which are spaces of meromorphic functions of degree N over Riemann surfaces of genus g) in 
genus g >2 (we skip the case g = 1 when the classification of stable bundles is rather special PP, see 
also remark 121 below). 

As before, denote the critical points of a meromorphic function vr(P) on a Riemann surface of 
genus S' by Pi,..., Pm and their images in CP^ by Im = '^{Pm)- This meromorphic function realizes 
the Riemann surface as a ramified A^-fold covering L of genus g of the Z-sphere with ramification points 
Pi,..., Pm- Assuming that all the branch points Im are different and simple, we get, according to the 
Riemann-Hurwitz formula: 

M = 2g + 2N -2 . 

Let’s introduce on L a canonical basis of cycles {aa,ba) (a = 1,...,^) and corresponding basis of 
holomorphic 1-forms Wa{P) {a = 1 ,... ,g), normalized by Wf^ = 6 ^/ 3 . 


3.1 Rank 1 systems in arbitrary genus 

The scalar systems (j2.4bll admit natural generalization to the Hurwitz spaces Pg, 7 v. Denote the prime- 
form on L by P(P, Q) (where P,Q € L); introduce the Bergmann kernel P(P, Q) = dpdg lnP(P, Q). 
By Bm{P) we denote the meromorphic differential of 2nd kind with vanishing a-periods and single 
pole at Pm of the second order with the following local behaviour: 

Bm{P) = f ^ + 0(1)) dxm 
\^m J 


as P ^ Pm, where Xm = \/Z — Im is a local parameter near Pm- Differentials Bm are related to the 
Bergmann kernel as follows: 


dXmiQ) 

Corresponding abelian integrals we denote by 0m- 


Q=Prr 


(3.1) 


0rr*(P) = 


Bm — 


d 


'Qo 


E{P,Q) 
dxmiQ) E{Qo,Q) 


In 


Q = Prr 


(3.2) 


where Qo £ T is a base-point such that its projection on Z-plane does not depend on {Z„}. 
Let us prove the following variational formula. 


Theorem 5 Assume that the local parameters xp, xq don’t depend on (some) Im- 
symmetric function 


b{P,Q) = 


B{P,Q) _ 52 


■lnP(P,Q) . 


dxpdxQ dxpdxQ 

Then dependence of b{P, Q) on the branch point Im is given by the following equation: 


Introduce the 
(3.3) 


db{P,Q) 

dim 


^b{P,Pm)b{Q,Pm) 


(3.4) 


16 











Proof. Formula (imi) is closely related to the Rauch variational formulas IZilEl which describe de¬ 
pendence of holomorphic differentials on the branch points. The proof is also very similar. Namely, 
consider the /^-derivative of the Bergmann kernel: dB{P,Q)/dim- This is a symmetric 1-form on 
L X L. Consider the Taylor series of B{P, Q) with respect to its hrst argument P in a neighbourhood 
of ramification point Pm- 

B{P, Q) = {ao -b aiXm + ■■■}dxm = {oo -b aiXm + ...}—, (3.5) 

zy L 

where oq, ui,... are some 1-forms with respect to Q. Differentiation of (13.51) with respect to Im gives 


dB{P,Q) 

dl 


{ao -b o(l)} 


dl 

4(/ - Imf/^ 


{ao -bo(l)} 


' 


(3.6) 


as P —> Pm- Therefore, dB{P, Q)/dlm is a meromorphic 1-form with respect to P, with the only pole 
at Pm of the second order with leading coefficient 


ao{Q) 

2 


B{Q,P) 


2dXm{P) P=P 


(which is itself a 1-form with respect to Q) and vanishing a-periods. Taking into account the symmetry 
of B{P, Q), we get 

dB{P,Q) lf P(P,P) 1 f B{Q,R) 
dim 2 ^ dX-mi^R) R=P7n J \ dXm(^) 

which is equivalent to (TOl) . □ 

Since differentials Bm are related to the Bergmann kernel via (HU), the theorem immediately 
implies the following useful 



Corollary 2 The abelian differential BffP) and the abelian integral 0n{P) depend on Im (for any 
m n) as follows: 


d 

dim 


Bn{P) 


-b{Pm, Pn)Bm{P) , 


^0n{P) = lb{Pm,Pn)0miP) 

Olm ^ 


(3.8) 


Now we are in position to formulate the analogs of Euler-Darboux equation in arbitrary genus: 


Theorem 6 Consider the following linear system of sealar equations for function 'ip{P, {/m})-' 

= Rm0m{P) , (3.9) 

atm 

where Rm are some funetions of {Im}- Denote by Pq a point of L such that Iq = vr(Po) does not depend 
on {Im}- Then the compatibility conditions of the system id.tA) are given by the following equations: 


d^f _ b{Pm,Pn) f ^ 

dlmdln 2 { Vm dim 'O-n dff j 


where f{{lm}) = ip{Po) and 


Om — 0m(,Po) 


funetion b{P,Q) is defined by 
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Proof. Substituting in (Eini) P = Pq, we see that 


Pm. — 


f 


0m{Po) ' 

Then the compatibility conditions of the system dSSI) are equivalent to the system of equations 


g r df 0n{P) ] d fdf 0miP) ] 
dlm\dln0n{Po)i Bln \ dim 0m{Po) j ’ ^ 

Using (Tm . we rewrite this equation as follows: 


(3.11) 


[ 0n{P) 0m{P) ] 

\0n{Po) 0m{Po)i 


{0„(P)0„(P„) - 0„(P)0„.(P„)) = 0 . (3.12) 

The l.h.s. of this equation is an abelian integral on L with respect to argument P with vanishing 
a-periods (since all o-periods of differentials Bm{P) vanish) and the poles (of the first order) only at 
Pm and Pn- Since 0m{P) = — 1/a^m + 0(1) as P —> Pm, we immediately see that equations (Tnni) 
are equivalent to absence of poles of at Pm and Pn- Therefore, the l.h.s. of (ItTI) must be a 

constant with respect to P; choosing P = Pq we see that this constant vanishes. We conclude that 
equations imni) indeed provide compatibility of the linear system (TUH) . 

□ 

In analogy to the case of genus zero, now we shall define the tau-function of the system (imH) and 
construct its solutions via solutions of the scalar Riemann-Hilbert problem on L. 

We define the tau-function of the system (IXTUl) by the following system of equations: 


A. 

din 


f2 

■Inr = -^ , 


din' 

dim. 


= 0 


From variational formulas (EH) and equations (imH) it follows that 



= b(Pm,Pn) 




VmVn 


(3.13) 


(3.14) 


Symmetry of this expression with respect to m and n proves compatibility of equations Eiini). 

As well as in the case of genus zero, we can prove the following 

Lemma 1 Definition iV. of the tau-function can be alternatively rewritten in the following form: 


cllnr 

dim 


1 

-res 


Proof. Choosing the standard local parameter Xm 
therefore, 


res 


I (dp^l 



dl J ■ 

= - Im near Pm, 



we have dxm 


(3.15) 
dl/^Xm] 
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On the other hand, as P belongs to a neighbourhood of P^, the linear system (ESI) can be rewritten 
as follows: 

dijj _1 <9V’ _ df 0m{P) 

dim dXm dim 0m(-Po) 

where we separated the dependence of V’ on Im which comes from /^-dependence of Xm- Taking the 
residue at Pm, we get 

d'>P X _ o flm 
dxm^ 0m(Po) ’ 

which implies coincidence of (Trni) and (jd.lSI) . 

□ 

The next theorem provides solutions of the system innni) . 


Theorem 7 Let I be an arbitrary smooth closed contour on L such that its projection on the I-plane 
7 r(/) is independent of {Im} and Pm 0 I for any m. Consider on I an arbitrary Holder-continuous 
function h{Q) independent of {Im}- 
Then the function 

satisfies the system UlM . Corresponding solution of the linear system /TO) is given by 

= (3.17) 


Proof. Before proving the statement of the theorem we notice that if the function h(Q) is Holder- 
continuous, then (UrTTl) is a solution of the Riemann-Hilbert problem on the contour / with jump 
2'Kih{Q). 

Let us now verify that the functions / and if defined by (ITTHl) and (IXT71) . satisfy the linear system 

(TnUl : 

I _ &miP) r /Q 1 

“ 0m{Po)^^^ ■ ^ ^ 

Taking into account representation dH of 0m{P) in terms of prime-forms, in analogy to Rauch 
formulas given by ThO we see that the Cauchy kernel dq In depends as follows on Im (assuming 

that all points P,Q,Qo are /^-independent): 

^ E{Qo% } = 2^^^Q)®rn{P) - (3.19) 


Therefore, 




= ^ j>^HQ)Bm{Q) ■ 


0m{P) 

Independence of this expression of P proves (ITTHl) . Therefore, the function / satisfies the compatibility 
conditions of (IXTHl) i.e. the equations (imH) . 


□ 


We notice that, using superposition principle for the systems we can substitute integration 

over contour in (IxDB by integration with respect to an arbitrary {/m}-independent measure with 
compact support; the only condition we need to impose on this measure is commutativity of integration 
in ll3T6l) with differentiation with respect to {Im}- 
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3.2 Systems of higher rank 

Here we consider coverings of genus g > 2. To define a natural rank r analogs of equations it 

is natural to start from generalization of the zero curvature representation (CTl) to higher genus: 


(9T 


dir. 




(3.20) 


where the r xr “Lax matrix” Um{P, {In}), T" G L has only one singularity on L, which is a simple pole 
at the ramification point Pm- We shall assume that trUm = 0 (this condition can be always satisfied 
by normalization of function T: T ^ [det'L]“^/”'I'). As well as in rank 1, the “Lax matrix” Um{P) 
can not be single-valued on L. In rank 1 matrix Um has additive twists along basic cycles on L, but 
dpUm{P) is a single-valued meromorphic differential on L with pole of second order at Pm- 

In higher rank this is not enough - it is necessary to introduce more degrees of freedom, and allow 
dpUm to suffer similarity transformation under tracing along each topologically non-trivial closed 
contour on L (this transformation may be {/m}-dependent itself). 

More precisely, consider a stable vector bundle y of rank r and degree d over L. Require dpUm to 
be a meromorphic section of the bundle ad^ ® K with quadratic pole at the ramification point Pm {K 
is the canonical line bundle), i.e. 

dpUm G H^{L, ady- (g) K{2Pm)) ; 

according to terminology proposed by Hitchin m, dpUm is called the meromorphic Higgs held. 

According to Narasimhan-Seshadri theorem m, the stable bundle y is characterized by the set 
of unitary constant (independent of a point P G L) matrices Xai , Xag, Xh, ■ ■ ■, Xbg, satisfying one 
relation 

d. 


n Xa,XbgXa-Xb. = exp{27ri-} . 


a=l 


Stability of the bundle y implies h^{ad^) = 0 i.e. the bundle ad^ does not have any holomorphic 
section. On the other hand, hP{ad^ (g) K) = (g — l)(r^ — 1); therefore, hxing the singular part of dpUm 
at the ramihcation point Pm'- 


R 


dpUmiP) = ( -^ + 0(1) ) dx 


as 


P^Prr 


(3.21) 


we dehne dpUm up to an arbitrary linear combination of {g — l)(r^ — 1 ) holomorphic sections of 
ady- ® K. In rank 1 case we hxed dpUm uniquely imposing the normalization condition of vanishing 
of all o-periods of dpUm', unfortunately, such simple normalization in higher rank is not known to the 
authors. 

The “Lax matrix” Um is uniquely dehned by dpUm inside of the fundamental polygon L of L if 
we hx a normalization point Pq G L and assume that Um{Po) = 0- The only singularity of the matrix 
Um is the simple pole at Pm of the form 

UmiP) = — + 0{l) as P^Pm', 

Xm 

generically Um has both multiplicative and additive non-single-valuedness along any cycle 7 G 7 ri(L): 


UmiP'^) = X^UmiP)Xy + Cl 


(3.22) 
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since dpUm{P'^) = XydpUmiP)X'y^ ^ the “additive twists” Cm don’t depend on P. 

The compatibility conditions of the linear system are given by: 

Fmn{P) = 0 , (3.23) 

where Fmn{P) is the curvature: 

Fmn{P) = ^-^ + [Um,Un]. (3.24) 

Ol'ji Oijyi 

To rewrite the compatibility conditions in terms of variables depending only on {Im} (and not 

on the point P of the covering L), we consider the coefficients of the Taylor series of Um{P) at the 
ramification point Pn, n ^ m: 


Um{,P) — Smn T TmnXn T 0(^n) • 

Then non-singularity of Fmn at Pm is equivalent to the condition 

+ ^Tnm + [Rm, Snm] = 0 , m^n . (3.25) 

The non-singularity of Fmn is insufficient for its vanishing, since Fmn{P) has both multiplicative and 
additive twists along topologically non-trivial loops. However, if we require that the “additive twists” 
Cm from (I3.22|) are related to matrices X') via equations 

= CZ, , (3.26) 

we observe that transformation of the Lax matrices Um is nothing but the gauge transformation 

of connection 1-form 'YZm=iP'rndlm by the matrix X'^ii^m})- Then the curvature coefficient Fmn 
transforms as follows: 

Fmn{P^)=X^Fmn{P)X^^ 

i.G. Fyfiji G H thus -^rn.n — 0* 

We conclude that the compatibility conditions of the linear system in arbitrary rank r and 

any genus g >2 are given by the system of equations and (TM . The new feature in comparison 

with the genus zero case is that we get another degree of freedom - the stable bundle x must itself 
depend on the branch points Im according to equations (since generic vector bundle over L is 

stable, generically evolution preserves stability of x at least locally, in a neighbourhood of a 

given stable bundle). 

Obviously, without any normalization of the twisted 1-form dpUm the coefficients Smm Tmn and 
additive twists Cm are not uniquely determined by the set of residues Rm and matrices X')- Therefore, 
the number of equations is substantially smaller than the number of variables. 

A possible way to define dpUm uniquely is to make use of one of meromorphic bidifferentials 
W{P, Q) on L X L whose existence is provided by the following lemma: 

Lemma 2 There exists meromorphie bidifferential W (P, Q) on L x L satisfying the following condi¬ 
tions: 
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1. On the diagonal P = Q the bidifferential W{P,Q) has second order pole with biresidue equal to 
X matrix 11 (which is the permutation matrix in C’’ <S> C^J: 


W{P,Q) = 


n 


{xp - XqY 


+ 0(1) > dxpdxQ . 


2. Symmetry condition: 


W{P,Q)=mV{Q,P)Ii . 

3. Automorphy conditions: for any 7 G vri(i) we have 


W{P"f,Q) =X^ W{P,Q) X. 


-1 


7 5 


(3.27) 

(3.28) 

(3.29) 


W{P,Q^)=X^W{P,Q)X^ , (3.30) 

1 2 

where for any linear operator A in PA we denote by A and A the operators A® I and I ® A in 
( 8 ) , respeetively. 


The relation (TOnil is obviously a corollary of and the symmetry requirement (TT^ . Equiv- 

alently, relations (IX^ . mean that W{P,Q) belongs to H^{ad^ ®K (2Q)) with respect to its 

2 

first argument, and to H^{ad^ (2P)) with respect to its second argument. 

Proof. Existence of bidifferential Wo{P, Q) satisfying only (13.2711 and automorphy properties (TM . 
(ICTl) . can be proved similarly to, say, existence of Schiffer kernel corresponding to the bundle ad^ (see 
for example 0 )- To construct bidifferential W(P, Q) which satisfies in addition the symmetry condition 

dSUHl), suppose that F{P, Q) = Wq{P, Q) —niTo(Q) T*)!! does not vanish; this is a holomorphic section 
1 2 

of ady^ <SiK with respect to P and holomorphic section of ad^ with respect to Q. Obviously, 

nF(p,g)n = -F(Q,p). 

Now define W{P,Q) = Wq{P,Q) — ^F{P,Q). Simple calculation shows that it satisfies (I3.28|) : other 
required properties are inherited from Wo{P,Q). 

□ 

We proved existence of W{P, Q); obviously this bidifferential is not unique: we can add to W{P, Q) 
an arbitrary linear combination of bilinear products of holomorphic sections fk{P) of ady- (8> K, satis¬ 
fying the symmetry condition dSUHl); this is a linear combination of the form 




E 

j,k=l 


Oijk 


fk (P) fi iQ)+ /, (P) fk (Q) 


with arbitrary Ojk G C. 

Let us fix the bidifferential W{P,Q) in some way (for example, according to remark |31 below). 
Then the Higgs fields dpUm can be defined as follows: 


dpUm{P) 


12 


tr W (P, Q) Rr 


dXm{Q) 


Q=Prr 


(3.31) 
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Now all variables Smn, Imm Cm become functionals of the residues Rm, branch points Im and 
matrices X'y] the number of variables {x^y^Rm) in the system (ld.25ll . (I,‘'{.2hjl coincides with the number 
of equations. 

The natural definition of tau-function, in agreement with (EUl), looks as follows; 


d 

dim 


Inr = 


1 

-res 


tr(dpTT“i)2 

dl 


Consistency of the definition is provided by the following 


(3.32) 


Lemma 3 Let function T solve the linear system with Lax matrices Um satisfying condition 

Um{Po) = 0 for some Pq G L, and conditions \d.dl\) . Then the 1-form 



tr(dp^'T“i)2 
Pm dl 


dim 


(3.33) 


is closed. 


Proof. Consider the Taylor series of T(P) in a neighbourhood of Pm- 

T(P) = ^0 + Xm^i + 0{xl^) ; 


then 

and 

Therefore, 

and 


1 I 

-res I 


'k.m'h-' = 'hid/o^ + 0(l), 




Rm = 


tr((ip'h'I' ^)^ 1 I tr{—2Rmdxm)^ 


dl 


2 res I 


2xmdx. 


= tiRz 


m^'^m 


Now we have to make sure that the derivative ^tri?^ is symmetric under the interchanging of m 


and n: 


A 

din 


As P —> Pmi we have: 


trP^ = -2iTRm { -Trim + [Rm-, Snm] \ = -tlRmTum - 


dXr, 


(3.34) 


tl {dpLfm{P)dpLfn{P)} = {-tlRmTnm + --- 

symmetry of (TOl) is thus equivalent to relation 

bires|p^tr{dpC/mdpf7n} = hires\pjT{dpUmdpUn} - 

Let us rewrite the l.h.s. of this relation in terms of bidifferential W{P, Q) = w{P, Q)dxpdxQ according 
to (ITTIT) . taking into account the behaviour H3.27I) of W{P, Q) on the diagonal P = Q: 


1(3,3 13 


2 , 2 12 


bires|p ti{dpUmdpUn} = bires|p tr tr { RmW {P,Pm)} tr { RnW {P,Pn)} ^ {dxm{P)) 
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Similarly, 


1 2 r 1 2 12 

—tr tr s RrnRn^ {,Pm^Pn) 


I 1 2 r 2 1 12 

birGs|p tl^dpUradpU'ii'^ — trtr s RmRn'^ i,Pn: Pm') 


(3.35) 


12 12 1 


12 


— tr tr < RraRn n W j Pji)Tl > — tr tr < RmRn'^ (Pmj Pfi) 


12(1 2 12 


(3.36) 


coinciding with (E3SI); the last equality in follows from the symmetry property (13.2811 of 

W{P,Q). □ 


Remark 2 More explicit treatment is, as usual, possible for elliptic coverings g = 1, when the bundle 
ad^ can possess meromorphic sections with single simple pole. In this case the additive twists Cm are 
absent, and monodromy matrices X 7 of the bundle x can be chosen to be independent of {Im}- la this 
case the above construction can be nicely rewritten in terms of elliptic r-matrix m- 


Remark 3 In rank 1 case we have fixed the bidifferential W{P, Q) by the requirement that all of its 
a-periods vanish; then W{P,Q) coincides with Bergmann kernel. This kind of normalization is not 
possible in higher rank due to non-invariance of W(P, Q) under tracing along topologically non-trivial 
loops on L. However, in rank 1 case there exists another way to fix W{P, Q) uniquely: one can require 
that v.p. ffj^W(P, Q)wa(Q) = 0 for any holomorphic differential Wa (in this case W(P,Q) is called 
the Schiffer kernel). These conditions have natural higher rank analog [SI. Namely, we can require that 
2 12 

v.p. ffj^ tr {W {P, Q) fl (Q)} = 0 for any G P[^{ad^ 0 K)] due to unitarity of the matrices Xa the 
integrand is a (1, l)-form on L. If the biresidue at P = Q is chosen to be the unit matrix 101 instead 
of permutation matrix H, our definition of W(P, Q) would give the Schiffer kernel corresponding to 
the bundle ad^. [S]. This normalization of W{P,Q) leads to the following normalization of dpUm- 
^■P-ffh = 0 for any G H^{ad^ 0 K). 

We notice that, in contrast to Bergmann kernel in rank 1, this normalization of W{P,Q) makes 
it non-holomorphic function of “moduli” lm\ therefore, in principle, the complete system of equations 
should contain also equations with respect to {Im} We shall discuss these aspects in more details in 
further publication. 


Relationship to isomonodromy deformations in higher genns. The link between the genus 
zero systems (j2.1311 and isomonodromic deformations on Riemann sphere discussed in Sect 12.51 can be 
extended to arbitrary genus (we outline it here for g >2). Let us briefly describe the isomonodromic 
deformations on the covering L. Consider a stable bundle x characterized by a set of unitary matrices 
X 7 for any 7 G vri(L); consider also a divisor Q = Qi + ■ ■ ■ + Ql on L. 

Introduce a Higgs field A{P) which is allowed to have simple poles at the points Qi, ... ,Ql, i-e. 
A G H^{ady- 0 K{Q))] suppose that trH(P) = 0. The higher genus analog of the linear differential 
equation looks as follows: 

= H(P)T . (3.37) 

where function T has unit determinant and satisfies the initial condition 'k(Po) = I at some point 
Po GL. 
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In analog to genus 0 case, function 'I' has regular singularities at the points Qi, ■ ■ ■ ,Ql with some 
monodromy matrices i.e. under tracing around Qk the function ^'(P) transforms as follows: 

^'(P) ^ ^(P)Mfc . 

Under tracing along basic cycles of L the function 'h gains left multipliers given by the matrices x^'i 
addition, it may gain the right multipliers Ma^, which are analogs of monodromy matrices M^: 

= Xa.^(P)M,„ , 'I/(P'’“) = y,„xh(P)M,„ . 

The monodromy matrices M^., Ma^, generate a SL{r) “monodromy” representation of the fun¬ 
damental group 7ri(L \ {Qfc})- 

If we now assume that all the monodromy matrices M^, Ma^, M(y^ are independent of the branch 
points Im and projections = '^{Qk) of the regular singularities, then function T satisfies the defor¬ 
mation equations 


= 144 -, 

( 3 . 38 ) 

= Um^ , 

( 3 . 39 ) 


where matrix 14 has simple pole at the regular singularity Qk\ matrix Um has simple pole at the 
ramification point P^; these matrices transform as follows under the tracing along any 7 G 7 ri(L); 

Pfe(p^) = x^Vk{P)x:;^ + (x7)/..x:7' , 

Um{P^) = X^Um{P)x:^^ + (X7)z^X7 ' • 

Obviously, the part of the deformation equations is nothing but the linear system (ItM in¬ 

troduced above; therefore, the isomonodromic deformations in higher genus correspond, as well as in 
genus zero, to a subset of solutions of the integrable systems (IT^ . (IT^ . 

4 Systems of rank 1, Darboux-Egoroff metrics and systems of hy¬ 
drodynamic type 

4.1 Darboux-Egoroff metrics 

It turns out that each solution of the rank 1 system (IXTUl) defines a flat diagonal (pseudo-)metric 
in (which gives rise to a flat diagonal Darboux-Egoroff metric in if additional reality and 
positivity conditions are imposed). In the sequel we shall (in agreement with previous works on the 
subject) use the term “Darboux-Egoroff metric” for such pseudo-metrics in C'^. 

For diagonal (pseudo-) metric 

M 

ds — ^ 9mm, dl^ ( 44 ) 

m=l 

the Christoffel symbols are given by 

^mn ~ 0 ) ^nm ~ \/ 9nn i m ^ Tl ^ k ] ( 4 - 2 ) 
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they are related as follows to rotation coefficients jSmn- 


„ 'J9mm V 9mm , 

Pmn = Tmn = - 7==— > m / n . 

Y ynn Y ynn 

The metric gU) is flat iff the rotation coefficients satisfy the following equations: 

dPr, 


dh 


— PmlPlr 


(4.3) 


(4.4) 


for any distinct l,m,n, and each of Pmn is invariant with respect to simultaneous shifts along all {Ik}- 


M 




k=l 


dlk 


(4.5) 


If in addition the rotation coefficients are symmetric, Pmn = Pnm, which is equivalent to relation 
draQuri = dngmm, then there exists potential U such that 


9mm — 


dU 

dim. 


and the metric (EU) is called the Darboux-Egoroff metric. 

We shall now prove that each solution of the system (Tnnii corresponds to Darboux-Egoroff metric. 


Theorem 8 Let f be an arbitrary solution of the system 123); and T{{lm}) be the corresponding 
tau-funetion defined by 1221 or urm) . Then metric with 

cllnr 

gmni — (^• 6 ) 

is a Darboux-Egoroff metric in . 


Proof. Let us compute the rotation coefficients of the metric EtJ. From the definition of tau-function 
(Trni) we have 


d\n.T 

dim 


Vm 


using the expression (IXTH) for the second derivative of the tau-function, we find 


fj _ _ (fl^ _ _ \h( P P 

■ 2 v/(h^y(I^ 


(4.7) 


These functions satisfy the equations dU as a corollary of variational formulas 

It remains to prove that each b{Pm,Pn) satisfies equations (14.511 . We shall use invariance of the 
Bergmann kernel B{P, Q) with respect to biholomorphic maps. Let us consider the branched covering 
L*^ which is obtained by small e-shift of all the ramification points Pm in /-plane i.e. the projections of 
branch points Pf^ of L*^ on /-plane are equal to /^ = /^ -|- e; is the Bergmann kernel on L*^. Denote 
the projections of points P and Q on /-plane by / and g, respectively. Define the point P^ to be the 
point lying on the same sheet as P and having projection / -|- e on /-plane; in the same way point Q'' 
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belongs to the same sheet as Q and has projection + e on /-plane. Since can be holomorphically 
mapped to L by transformation / —> / -|- e on all the sheets, we have 

B\P^,Q^) = B{P,Q). (4.8) 

Assuming that P belongs to a neighbourhood of the branch point P^-, and Q belongs to a neighbour¬ 
hood of the branch point we can write down the respective local parameters as Xm{P) = — Im 

and Xn{Q) = y/— In- These parameters are obviously invariant with respect to simultaneous e-shifts 
of all {/m}) I and /j.: = Xm{P) and x^{Q^) = Xn(Q)- Therefore, equality (jd.SI) induces the 

same relation between b{P,Q) = B{P,Q)/dx{P)dx{Q): 

P{P^,Q^) = biP,Q) . (4.9) 

Assuming now that P = Pm and Q = Pn and differentiating (lOll with respect to e at e = 0, we come 
to H4.5f) . 

□ 

Of course, besides simultaneous translations, the Bergmann kernel is invariant with respect to 
any other Mobius transformation of /-plane performed simultaneously on all the sheets of L. We can 
use this invariance to obtain two relations, corresponding to other one-parametric families of Mobius 
transformations. 


Proposition 1 

ing relations: 


The rotation coefficients of the Darboux-Egoroff metric satisfy the follow- 

( M 


T,k 

s k=l 


A. 

dlk 


Pmn — /^r? 


(4.10) 



— ~(/n + / m )Pn 


(4.11) 


Proof. Relation (ICT71) corresponds to invariance of the Bergmann kernel under simultaneous dilatation 
on every sheet of L i.e. to the transformations 


^ (1 + e)^m ; / —> (1 -|- e)/ , /i —> (1 -|- e)fi . 

The new feature in comparison with the proof of relation (1131) is that the local parameters are now 
dependent on e: 

Xm{P^) = [(1 -h e)l - (1 -h = (1 + A^‘^Xm{P) 

and 

= [(1 + - (1 + = (1 + ; 

therefore, invariance B^{P^,Q^) = B(P,Q) of the Bergmann kernel translates on the level of b{P,Q) 
as follows: 

{l + e)P{P^,Q^) = b{P,Q) . (4.12) 

Differentiating this relation with respect to e at e=0 via the chain rule and choosing P = Pm and 
Q = Pn, we get (I4.10() . 
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In a similar way we can deduce (BTTTl) from invariance of the Bergmann kernel with respect to the 
one-parametric family of transformations 


I 




1 + el 


on each sheet of L. □ 

Relation (1071) for the rotation coefficients can be found in [HI; relation dnu) seems to be new. 
We also notice that all primary differentials used in to construct Frobenius manifolds from Hurwitz 
spaces can be obtained from solutions (lOl) of the systems (1071) by appropriate specification of the 
contour I and the function h{Q). 


4.2 Systems of hydrodynamic type 

According to well-known results (see, for example, review mi to each Darboux-Egoroff metric one 
can associate a class of diagonal systems of hydrodynamic type 


dx dt 


(4.13) 


where functions Vm{{lk}) (the “characteristic speeds”) are related to Christoffel symbols (14.21) of the 
metric 631 via system of differential equations: 


dmVn 


KmiVm - Vn) . 


(4.14) 


Compatibility of equations dm is provided by equations (lOl) for rotation coefficients. 

Let us choose some solution of the system (fTTITl) parametrized by an arbitrary function h{P) on 
contour I dm. Then the metric coefficients dm are given by: 

ffmm = fh{Q)Bm{Q)^ ■ (4.15) 

The Christoffel symbols of this metric look as follows: 


, §ih{P)Bm{P) 


(4.16) 


Solutions of equations for these Christoffel symbols are described by the following proposi¬ 

tion^. 


Proposition 2 Let hi{P) he an arbitrary Holder-continuous and independent of {Im} function on 
contour 1. Then the functions 

_ §,h{P)Bm{P) 

™ §,h{P)Bm{P) 

satisfy system with Christoffel coefficients given by O- 

^Formula lITTTIi is due to T.Grava (private communication); similar formula for the case of Whitham equations was 
obtained in EHl 
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The proof of this proposition is a simple calculation based on Rauch variational formulas for differen¬ 
tials Bm'- 

= l3mnBn{P) (4.18) 

Otn 

To construct solutions of the system with characteristic speeds gizi one needs to use the 

following theorem by Tsarev m-- 


Theorem 9 Let functions ■ ■ ■, Im) satisfy equations 


Then system of equations 


$m({4}) = i + Vm{{lk})x 


(4.19) 


defines implieit solution {lmix,t)} of the system of hydrodynamic type where ■ ■ ■, Im) is 

an arbitrary solution of the system of differential equations 






dnVm 

Vm-Vn 


(4.20) 


for m,n = 1,..., M. 


To apply the hodograph method to any of these systems we need to solve also the system of 
equations for functions Obviously, the solution is given by same formulas as solution 

of equations (1331 for functions Vm, but with another arbitrary Holder-continuous function h 2 {P) 
independent of {Im}'- 


§^h2{P)Bm{P) 

§,h{P)Brn{P) • 


(4.21) 


For each choice of h 2 {P) the system of equations (14.191) defines the implicit solution {lm{x,t)} of the 
system of hydrodynamic type (jUl- 
We proved the following 


Theorem 10 Consider system of hydrodynamic type where velocities Vm are given by formula 

with arbitrary Holder-continuous and independent of {Im} functions h{P) and hi{P) on contour 
1. Let h 2 {P) be another arbitrary and independent of {Im} Holder-continuous function on contour 1. 
Then system of M equations for M variables {lm{x,t)}^^i 


j){h 2 {P) + h{P)thi{P)x}Bm{P) = 0 , m = (4.22) 

defines implicit solution {lmix,t)} of the system 

As before, the condition of Holder-continuity of functions h, hi and /i 2 on contour I can be relaxed. 
Namely, we can substitute the contour I by an arbitrary subset of L, and define on this subset three 
arbitrary measures independent of {Im}] the only requirement one needs to impose is commutativity 
of integration in with differentiation with respect to {Im}- 
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5 Summary and outlook 


In this paper we propose a new class of integrable systems of partial differential equations associated 
to spaces of generic rational maps of fixed degree. For maps of degree two such systems give rise to 
the Ernst equation from general relativity; for maps of higher degree our systems realize the scheme 
of deformation of autonomous integrable systems proposed by Burtsev, Mikhailov and Zakharov. We 
introduce the notion of the tau-function of the new systems, and describe their relationship to the 
matrix Riemann-Hilbert problem and the Schlesinger system. 

We generalize our construction to derive integrable systems associated to arbitrary Hurwitz spaces 
of meromorphic functions of degree N on Riemann surfaces of genus g >2. 

When the matrix dimension equals 1, our systems are linear; they can be solved via scalar Riemann- 
Hilbert problem on the Riemann surfaces. Each solution of such system corresponds to a flat diagonal 
metric (Darboux-Egoroff metric), together with corresponding systems of hydrodynamic type and 
their solutions. 

Our results suggest the following directions of future work. We expect our systems for g > 2 
to be natural deformations of two-dimensional version of Hitchin systems proposed in |2()j . The 
isomonodromic deformations in higher genus briefly discussed here should be in close relation to 
existing frameworks as well as to non-autonomous Hitchin systems HU. All these links 

should be clarified. The applications of the new systems should also be studied, especially from the 
point of view of their potential relationship with structures of Frobenius type. 
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